A simple electrostatic model applicable to biomolecular recognition 



T. P. DoeriE and Yi-Kuo YiQ 

National Center for Biotechnology Information, National Library of Medicine, 
National Institutes of Health, 8600 Rockville Pike MSC 6075, Bethesda, MD 20894-6075 

An exact, analytic solution for a simple electrostatic model applicable to biomolecular recognition 
is presented. In the model, a layer of high dielectric constant material (representative of the solvent, 
water) whose thickness may vary separates two regions of low dielectric constant material (represen- 
tative of proteins, DNA, RNA, or similar materials), in each of which is embedded a point charge. 
For identical charges, the presence of the screening layer always lowers the energy compared to the 
case of point charges in an infinite medium of low dielectric constant. Somewhat surprisingly, the 
presence of a sufficiently thick screening layer also lowers the energy compared to the case of point 
charges in an infinite medium of high dielectric constant. For charges of opposite sign, the screening 
layer always lowers the energy compared to the case of point charges in an infinite medium of either 
high or low dielectric constant. The behavior of the energy leads to a substantially increased repul- 
sive force between charges of the same sign. The repulsive force between charges of opposite signs 
is weaker than in an infinite medium of low dielectric constant material but stronger than in an 
infinite medium of high dielectric constant material. The presence of this behavior, which we name 
asymmetric screening, in the simple system presented here confirms the generality of the behavior 
that was established in a more complicated system of an arbitrary number of charged dielectric 
spheres in an infinite solvent. 

PACS numbers: 41.20.Cv,87.10.Ca 



I. INTRODUCTION 



The proper functioning of biomolecular systems depends upon the aggregation of multiple molecules embedded in 
a high dielectric constant solvent (water). From the medical point of view, there are both normal complexes (such 
as ribosomes) and abnormal complexes (such as amyloid formations). Understanding the microscopic mechanisms 
involved in the aggregation process would illuminate both normal and abnormal states, and could aid the modification 
of existing complexes or the design of new ones. This work examines the electrostatic interaction, among the most 
important interactions in biomolecular systems. 

In previous research that developed a scheme for computing to known precision the energy and forces in a system 
of an arbitrary number of charged dielectric spheres embedded in an infinite solvent Q, an effect that was called 
asymmetric screening was observed. Namely, the magnitude of attractive electrostatic interactions was decreased 
(relative to point charges in an infinite solvent) while the magnitude of repulsive electrostatic interactions was increased 
(again, relative to point charges in an infinite solvent). It was speculated that this effect might aid biomolecules such 
as proteins in the adoption of correct conformations and in intermolecular recognition. 

This paper presents further studies of this effect in a simplified system that is amenable to complete and thorough 
analytic examination. The simplicity of the model is an advantage in this case because one wishes to examine in more 
detail an effect that is already known to occur in the more general and less symmetric system of spheres mentioned 
above. The system studied here can be considered a simplified model of two molecular surfaces during the process 
of binding or aggregation. Instead of spheres, consider two half-spaces, each with a single point charge embedded, 
separated by an infinite slab of high dielectric constant material (water, for example). If the dielectric constants are 
swapped, then one would have a model of, for example, a membrane in water. Separation of variables is used to 
obatin the potential, and from that the energy and the force between the two half-spaces. It is more convenient to 
use the surface charge method [7(j-[9[ to obtain the density of surface charge induced on the two surfaces. 
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FIG. 1: The most general situation under consideration. The shaded region is infinite in the x and y directions, has thickness 
2d in the z direction, and is filled with a material with dielectric constant eo- The origin is chosen so that the distance from the 
origin to each surface of the shaded region is d. The unshaded region entirely in the z > half-space is filled with a material 
with dielectric constant e\ and contains a charge q\ on the positive z axis a distance d + si from the origin and a fixed distance 
Si from the surface of the shaded region. The unshaded region entirely in the z < half-space is filled with a material with 
dielectric constant £2 and contains a charge q-z on the negative z axis a distance d + S2 from the origin and a fixed distance S2 
from the surface of the shaded region. 



II. THE GENERAL SITUATION 



Consider a slab of material of thickness 2d, infinite in the other directions, with dielectric constant eo sandwiched 
between two half-spaces filled with materials of dielectric constant e\ and £2 respectively. A charge q\ lies within the 
external material with dielctric constant e\ a distance s\ from the internal material (dielectric constant eo); a charge 
52 lies within the other external material (with dielctric constant £2) a distance S2 from the internal material and a 
distance Si + S2 + 2d from the charge q\. Place the origin of coordinates half way between the two charges. Place the 
z axis through the line joining the two charges, perpendicular to the surfaces of the internal slab of material, and with 
the positive z axis passing through the charge q\, as in Fig. [TJ Because of the symmetry of the system, cylindrical 
coordinates (p, 4>, and z) will be used. 

We wish to find the electric potential ($), the electrostatic energy (U), and the force (F) required to pull the 



3 



external materials apart. We begin by determining the potential in the general case. Azimuthal symmetry implies 
that the potential $ is independent of 4>. The symbols $ : an d $2 will be used to indicate the potential in the 
interior material, in the material entirely in the positive z region, and in the material entirely in the negative z region 
respectively. The boundary conditions are 

1. <f> — > as z — > ±00 

2. <J> (z = d)= = d) 

3. <5> 2 {z = -d) = <P {z = -d) 
4- £q— \ z=d - £1— \ z=d 

K f „»2 I — Fn^a I 



The appropriate general solution of Laplace's equation is 

00 />0O poo 

/ J m (kp)(ae kz +be~ kz )(csmm<p + dcosm<j))dk — > / J Q (kp)(ae kz + be~ kz ) dfc, 



because of the azimuthal symmetry. The appropriate form of the potential of a point charge at p = and z = z' is 

e -*l»-*'lj (fcp) dfc. 



The potential in the positive z region of exterior material is a solution of Laplace's equation plus the potential of 
the screened point source: 



f°° n f°° 

$1 = / B 1 (k)e- kz J {kp)dk + ^- e- k ^ z - d - Sl ^Jo(kp)dk, 
Jo £ i Jo 



(1) 



where boundary condition 1 has deleted one of the exponentials in the solution of Laplace's equation. Similarly, the 
potential in the negative z region of exterior material is 



e^ z+d+s ^J (kp)dk. (2) 



p CO p C 

$2 = / A 2 {k)e kz J (kp)dk + ^ 

Jo £ 2 Jo 

The potential in the interior material is 

/>OC 

$0=/ (A (k)e kz +B (k)e- kz )Jo(kp)dk. (3) 



Boundary conditions 2-5 determine the coefficients: 

B (k) = c k ^-^ ^ (£ ° + £l)(£ ° Z £ ^ ~ eM4d+S2) (gQ ~ fiXgo + £2)gi + 4e fc ( 2rf +^) £o£l g 2 . . 

-(£o-£i)ei(£o-£2)+e 4fcd ei(£ + ei)(eo + e2) 

An(k) 2c k <*-« — ) ^^^^ (£ ° + £2)qi + e " Sl (£ ° = £l)g2 (4b) 

° l J " -(£o-£i)(£o-£2) + e 4fc ' i (£o + £l )( £ o + £2) 1 j 

° W -(£o-£i)(£o-£ 2 ) + e^( £o + £l )(£o + £2) 1 ' 



A (k) = e*^-'!-^) 4efc(2d+S2)g o £ 2gi ~ e fc ( 4d +^)( £o + giXgp - £2 )g 2 + e fcsi (go - £l )(e + e 2 )g 2 , ^ 

2 ' -(£o-£i)£2(£o-£2)+e 4fcd £2 (£o+£i)(£o + £2) 



Not surprisingly, interchanging the indices 1 and 2 in the expression for B\ turns it into A 2 . 
The distribution of free charge (the two point charges) and the potential determine the energy: 



U = \J p/* = | *iO» = o,* = d + « 1 ) + ^(p = o, 



z = -d- s 2 ), (5) 
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where the primes on the potentials indicate that the potential of the point charge in the corresponding region has 
been subtracted out in order to avoid infinite self-energies. Substitution of Eq. (JTJ) , Eq. @, and Eq. ^ into Eq. (|5|) 
yields 

(e + £i)(£o + £ 2 ) Jo l-a ia2 e- ikd 2e x J 1 - ai a 2 e~ ikd 

.92/ e \ e oii-a 2 ) ,. , s 

2e 2 Jo 1 - a x a 2 e 4kd 

where a x = (e - £i)/(£o + £i) and a 2 = (e - s 2 )/{e + e 2 ). 

Because we imagine this situation to be a simplified model of two molecular surfaces separated by a layer of water, 
the force should be obtained by imagining that the charges are fixed with respect to the materials in which they are 
embedded, but the thickness of the interior slab is allowed to vary. In other words, the force we are considering is the 
negative of the derivative of the energy with respect to 2d: 

or in scalar form for the magnitude 

p = _ldU 
2 dd ■ 

Clearly, this simple model neglects any internal rearrangement of the molecules during the process of interaction, an 
effect that is believed to be important in many cases. However, while a model designed to capture the behavior of 
specific molecules would need to include such an effect, our purpose is only to investigate one particular interaction, 
the very important electrostatic interaction, and so this point is not a concern here. The force is 

p = 4 gig2£o f°° c - k(2d+Sl+S2 ) k l + aia 2 e- 4kd ^ 

(e +£i)(£o+£2) Jo (1 - ai a 2 e- ikd ) 2 



q2 roo e -k( 2si +Ad) k 

+ e7 a2(1 ~ ai) / (1 - ai a 2 e-^ d r 

oo e _fc(2s 2 +4d)£, 

£ 2 ~^ Jo (1 - ai a 2 e- 4kd ) 2 



-«i(!-«2)/ n z Z .-4^2 dfc - (7) 



We now examine two particular cases. 

III. TWO IDENTICAL CHARGES IN IDENTICAL MEDIA 

Let qi = q 2 = q, e\ = e 2 = e c , £o = £;, and si = s 2 = s. We are now considering a slab of material (thickness 
2d and infinite in the other directions) with dielectric constant e\ sandwiched between two half-spaces filled with a 
material of dielctric constant £ c . (Internal material is indicated by the subscript 'i', and external material is indicated 
by subscript 'e'.) A charge q lies in the external material a distance s from the internal material. An identical charge 
q lies in the other semi-infinite external material a distance s from the internal material and a distance 2s + 2d from 
the other charge. See Fig. [21 with the positive charge chosen. The potential, the energy, and the force follow upon 
making the appropriate substitutions in Eqs. |T][3]), Eq. (J6)), and Eq. (0 respectively. (Alternatively, it is a simple 
matter to set up and solve the boundary value problem for this particular situation.) 

Making the appropriate substitutions in Eq. ©, letting a = (e^ — £ c )/(£i + £ c ), and using the identity (4£j£ c /(£i + 
£ c ) 2 ) = 1 — a 2 , one finds the energy: 

U = t r e -*k. e- 2fcd (4£ i£c /(£ i - £c ) 2 ) + q(e- 4fcrf - 1) dfc 



£ c JQ 

^2 /-oo -, „,„2kd 



T_ I c ~2k(s+d) !-«e 
£ Jo 1 ~ ae 



e ^^z^dfc. (8) 
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FIG. 2: A simplified situation considered in detail. The charges are now of equal magnitute and are constrained to be the same 
distance from the origin. The cases of identical charges and of opposite charges are both considered. Both unshaded regions 
have the same dielectric constant, referred to as e e . The dielectric constant of the shaded slab is now referred to as e%. 



One may evaluate the integral by expanding the denominator in a series: 



u = q - 



J n— n 



n=0 

2 00 



E 



2e c ^ \ s + (n + l)d s + nd 



E 



2 



2e c a ^— ' s + nd 2e c as 

n— 

<7 2 (l-a 2 ) /s s \ g 



2e c o;s 



2^(1 l;£ + l;a)--^-. (9) 

\d d ) 2e e as 



where 2-F1 is a Gauss hypergeometric function. 

Even though the series in Eq. © was obtained by separation of variables, it can be interpreted as the effect of an 
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FIG. 3: Graphs of the energy as a function of separation, both for identical charges and for opposite charges. For comparison, 
the energy of point charges, both identical and opposite, in an infinite uniform medium (both e a and e{) is shown. The 
calculations are for e c = 1, ei = 80, s = 1, and q = 1. For opposite charges separated by a high dielectric layer, the energy 
varies little. For like charges separated by a high dielectric layer, the energy at small separations changes rapidly. 



infinite sequence of image charges. The charges have separations 2s + 2nd for n = 0, 1, 2, The magnitude of the 

image charges can be read off from the coefficients of q/(e c (2s + 2nd)) with appropriate care taken to separate out 
the direct interaction of the free charges. This interpretation brings to mind recent work that used an approximate 
series of image charges to study a pair of membranes in a solvent of water and ions. (Tlj 

Because the dielectric constant of water (« 80[I3|) is much larger than the dielectric constant of protein (w 4 13]), 
we are most interested in screening situation: < a < 1. In the limit a — ► 1, the interior slab becomes metallic. In 
this case we find that U = —q 2 /(e c 2s), which is just the interaction energy of each free charge with its image charge 
due to the metal; the two free charges do not 'feel' each other. If the media all have the same dielectric constant, then 
a = and U — q 2 /(e c 2(s + d)), which is simply the energy of two charges in an infinite dielectric medium. Similarly, 
if d = we find the obvious result U = q 2 /(e c 2s). Finally, in the limit that d — > oo, U — > — (q 2 a)/(e e 2s) < 0. In 
this case, the two fixed charges do not see each other, but each point charge can still induce a charge density on the 
nearby surface, and this process will always reduce the energy. Therefore U is negative in this limit. The behavior 
just summarized can be seen in Fig. [3] and Fig.HJ 



Making the appropriate substitutions in Eq. ([7]) and again using the identity (4£i£ e /(£i + £c) 2 ) = 1 



the force: 



F = 



{I -a 2 ) [°° fc e -2fc(d+s) 



(1 - ae~ 2kd ) 2 



dfc. 



one finds 



(10) 



Rather than performing a similar procedure with series to evaluate the integral, one may simply differentiate the 
series for U: 



F = — 



q 2 a(l-a- 2 ) 



q 2 (l-a 2 ) 



E 



H (2s + n2dy 

71 — V / 



4e c a ^— i (s + nd) 

n— ' 



2 ' 



(ii) 



As noted above, for the case of complete screening (i.e., a = 1) the free charges do not 'feel' each other. As expected, 
the force vanishes in this case. If the media all have the same dielectric constant, then a = and F = q 2 / (e c (2s+2d) 2 ), 
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FIG. 4: Graphs of the energy as a function of Ei, both for identical charges and for opposite charges. For comparison, the 
energy of point charges, both identical and opposite, in an infinite uniform medium (both e e and ei) is shown. The calculations 
are for 2s + 2d — 5, e c = 1, s = 1, and q = 1. 



the force between two identical charges in an infinite dielectric medium. On the other hand, if d = we find the 
curious result F = (g 2 £i)/(e 2 4s 2 ). When d = one might expect F not to depend on ej. However, F(d) samples U(d) 
in the vicinity of d, and even when d = a dependence is generated on £j, which characterizes the material that would 
fill the gap if one were to draw the two outer regions apart. Indeed, for d = and E\ — > 1, the force becomes infinite, 
i.e., the energy changes discontinuously at d = if a = 1. The behavior just summarized can be seen in Fig. [5] and 
Fig. © 

The difference between U and the energy of two point charges in an infinite medium of dielectric constant e c is 
defined to be AU. (This could not be calculated in the general case because in that case there is no single exterior 
material.) One finds 

„2 poo / 1 _ n,p 2kd \ n 2 n C°° 1 — p -ikd 

AU = — e~^ s+d ~> ( } a \ lr , -l)dk = - q — e- 2fc * i 6 Hfc. (12) 

Notice that AU < in the case of screening (a > 0), which makes sense because the energy should be lowered by 
replacing a portion of the low dielectric constant material with higher dielectric constant material. If a = 0, the 
energy U is the same as the term we have just subtracted off, so AU = 0. Similarly, if d = 0, then AU = 0. 

The force difference AF corresponding to AU can be obtained either from the expression for AU or the expression 
for F: 

AF = £^ te -« + .»(_il__^_ (13) 

In the case of d — we find that AF = 2g J(i_ a ) ■ If a = 1, then AF = —q 2 /(e c (2s + 2d) 2 ) < which, as expected, 
is just the term we subtracted off to form AF. Clearly AF = if a — 0. The behavior of AF for small but non-zero 
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FIG. 5: Graphs of the force as a function of separation, both for identical charges and for opposite charges. For comparison, 
the force between point charges, both identical and opposite, in an infinite uniform medium (both e e and £i) is shown. The 
calculations are for e c = 1, Si = 80, s = 1, and q = 1. The inset is a close-up of the three curves near the x axis for small 
separations. 
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FIG. 6: Graphs of the force as a function of e\, both for identical charges and for opposite charges. For comparison, the force 
between point charges, both identical and opposite, in an infinite uniform medium (both e c and £i) is shown. The calculations 
are for 2s + 2d = 5 e e = 1, s = 1, and q — 1. 
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a may be deduced from the series expression for F: 



AF 



> 



4e c 

9 

g 

4e c 



n(l — a 2 )a n 
^— ' a(s + n<i) 2 

n— 1 v 7 



E 



E 



n(l 



4e c (s + d) 2 

^,2 



^ a(ns + neQ 2 4e c (s + <i) 2 
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(1 



4e c (s + df 
4e e (s + tf) 2 2 



E 



a 



4e e (s + rf) 2 



2 ^ + ^ 2 )- 



When > 0, the repulsion between identical charges is stronger than the case when both identical charges are in 
one uniform medium with dielectric constant e c . Upon letting e c — > 1 (see Fig. [5] and Fig. [5]), we see that one can 
have a repulsion larger than in vacuum, a counter-intuitive conclusion. The origin of this behavior can be deduced 
by returning to Eq. ([6|), the energy for the more general situation first described. Setting q\ = 0, q 2 — q, and S2 = s 
but retaining distinct dielectric constants in each region, we find 



U = 



2e 2 



e- 2ks (e- 4kd ai 



a 2 ) 



aia2e 

Oil 



4kil 



dk 



Q'2 



2s + 4(n+l)rf 2s + And 



and 



2 00 

^ = rE«> 2 

n—0 



(n + l)ai 



nu2 



(2s + 4(n + l)d) 2 (2s + And) 2 



Each factor of a\ (a 2 ) indicates an image reflection across the surface of the material with dielectric constant e\ (e 2 ). 
Notice that the induced charge of the leading term (proportional to a.\) is the same sign as the free charge because 
the image charge is located on the low dielectric side of the interface. If the image charge were located on the high 
dielectric side of the interface (eo < E\ and £o < £2) then the induced charge would have the opposite sign leading to 
an attractive force similar to the more familiar case of a charge near a conductor. 

Now consider the energy and force differences (A£7 and AF) when the comparison is made to the interaction with 
the £j material everywhere. The energy difference in this case, AU, is 



AU 



-2k(s+d) 



ae 



2A-</ 



-2kd 



Q 2 f°° r -2Ms+d) ( g i ~ £ e) + a(g e e' 
EiS t 



dk 



-2kd 



e K e 2kd ) 



1 



ae" 



-2kd 



dk. 



(14) 



In order to understand the behavior of AU, we observe that AU(d = 0) = (g 2 (ei — £ c ))/(£i£ e 2s) > with equality 
when E\ = e c (i.e., a — 0). However, as d — > 00, AU —> — (q 2 a) /(2e c s) < 0. Evidently, for any positive a, AU is 
positive for small d and becomes negative for sufficiently large d. This behavior can be inferred from Fig. [3] Given 

that (£i/e c )(l - a) = (1 + a), AF is 



AF 



> 



9 



d) 



ke -2k(s+ 



ke -2k(s+d) 



(£i/e B )(l 



(1 - ae~ 2kd ) 2 

(l + «) 2 
(1 -ae- 2kd ) 2 



1 



dk 



dk 



(15) 



fc e - 2fc(;5+d) [(l + cv) 2 -l]dfc>0, 



which guarantees that AF > 0, as would be expected based upon Figs. O and O 
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IV. TWO OPPOSITE CHARGES IN IDENTICAL MEDIA 



Consider the same situation as in the previous section except that the two charges are of opposite sign. Namely, 
let 51 = q, q 2 = —q, £\ — e 2 = s B , £0 = £ b an d s i = s 2 = s. See Fig. ® with the negative charge chosen. The 
potential, the energy, and the force follow upon making the appropriate substitutions in Eqs. Eq. and 

Eq. ([7]) respectively. (Alternatively, it is a simply matter to set up and solve the boundary value problem for this 
particular situation.) 

Making the appropriate substitutions in Eq. ([6|), letting a = (ej — s e )/(si + e c ), and using the identity (4ej£ e /(£j — 
e c ) 2 ) = 1 — a 2 , one finds the energy: 



U 



-2kd 



(1 



a{e 



-4kd 



-2k(s+d) _ 



ae 



2kd 



1 



-2kd 



dk. 



1) 



dk 



(16) 



Again, we are most interested in screening situation: < a < 1. When a = 1 (perfect screening), we find that 
U = —q 2 /(e c 2s), which is just the interaction energy of each free charge with its image charge due to the metal; 
the two free charges do not 'feel' each other. If the media all have the same dielectric constant, then a = and 
U = —q 2 /(s c 2(s + d)), which is simply the energy of two charges in an infinite dielectric medium. Similarly, if d = 
we find the obvious result U = —q 2 /(e c 2s). Finally, in the limit that d — > 00, U — * —(q 2 a)/(e e 2s) < 0. In this case, 
the two fixed charges do not see each other, but each point charge can still induce a charge density on the nearby 
surface, and this process will always reduce the energy. Note that if a is close to unity (e.g., a water solvent), U varies 
little as d goes from to 00. The behavior just summarized can be seen in Fig.[3]and Fig. 2] Comparing Eq. (fTo| with 
Eq. ||5J), one sees that the series for U is the series for identical charges with an overall minus sign and the substitution 
a — > —a. 



Making the appropriate substitutions in Eq. and again using the identity (4ei£ c /(ei — £ c ) 2 ) 
the force: 



F = 



<7 2 (1 



fc e -2k(d+s) 

,-2kd\2 



(1 



dfc. 



1 — a 2 , one finds 



(17) 



As noted above, for the case of complete screening (i.e., a = 1) the free charges do not 'feel' each other. As 
expected, the force vanishes in this case. If the media all have the same dielectric constant, then a = and 
F = — q 2 /(e c (2s + 2d) 2 ), the force between two opposite charges in an infinite dielectric medium. If d = we 
find the somewhat non-obvious result F = —q 2 /(eiAs 2 ), the explanation for which is the same as in the case of 
identical charges. The behavior of the force in the case of opposite charges is more consistent with naive intuition: 
the force with a high dielectric layer is somewhere in between the force with low dielectric everywhere and the force 
with high dielectric everywhere. The behavior just summarized can be seen in Fig. [5] and Fig. [S] Comparing Eq. (JTTJ) 
with Eq. IjlOp. one sees that the series for F is the series for identical charges with an overall minus sign and the 
substitution a — > —a. 

The energy difference AU is now calculated along the lines used in the case of identical charges: 



All = — e - 2k ^ 
£c Jo 

2q 2 a rCO 



1 - 



1 



J2kd 1 



1 + ae 



-2kd 



dk 



(18) 



2k( s +d) sinh2fcrf dfc<a 







ae 



-2kd 



Since U(d = 0) = —q 2 /(e c (2s + 2d)) and U — > — q 2 a/(2e c s) as d — > 00, it is clear that AU should be negative (see 



Fig. [3|. As expected, the energy difference AU vanishes both for d = and for a 
interacts with its image charge, and therefore AU = —(q 2 d)/(2e c s(s + d)). 
Now consider AF for opposite charges: 



0. For a = 1, each charge 



AF = ^ ke- 2k{s+d ^ 

£ o Jo 



1 



(1 



,-2kd\2 



1 



dfc > 0. 



(19) 



The magnitude of the attractive force between opposite charges with a screening layer is always less than when both 
charges are in one uniform dielectric medium with dielectric constant e c . This agrees with intuition upon letting 
e c — * 1. As expected, AF vanishes if a = 0. Also, AF = q 2 /(e c (2s + 2d) 2 ) if a = 1, which confirms that there is no 
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force between charges that have a metal between them. For d = 0, the force difference AF = (q 2 a)/ (2e e s 2 (l + a)) 
depends on a for the reason noted in the case of identical charges. 

The energy difference when the comparison is made to the interaction with the ej material everywhere is AU: 



AU = 



3 -2/c(s+d) 



1) 



1 



Jlkd 1 



~2kd 



-2k(s+d) 



1 + ae 
(e c - £i) + a(e c e 



dk 

-2kd 



e K e 2kd ) 



1 + ae 



-2kd 



dk. 



(20) 



As expected on the basis of Fig. [3]), AU is less than or equal to since both terms within the square brackets are 
less than or equal to in the case of screening (0 < a < 1). For a = 0, the energy difference AU vanishes, while for 
a = 1 and d — > oo, AU = —q 2 /(e c 2s), the energy of interaction due to the presence of image charges. For d = 0, 
AU = q 2 (e c - £i)j2e e EiS. 
Now consider AF: 



AF = -2- 

£i Jo 

£ i Jo 



ke -2k(s+d) 
ke -2k(s+d) 



gj(l - q)(l + «) 
e e (l + ae~ 2kd ) 2 

(l + « ) 2 
(1 + ae 



- 1 



dk 



-2kd\2 



- 1 



dk < 0. 



(21) 



The attraction between unlike charges in our setting is always stronger than when the charges are in a uniform 
dielectric medium of dielectric constant e{. Clearly, AF vanishes when a = and when d = 0. 

V. COMMENTS 



The energy and force for the case of two point charges in a dielectric medium with a layer of differing dielectric 
between them has been compared with two baselines: point charges in a uniform medium having the dielectric constant 
of the separating layer and point charges in a uniform medium having the dielectric constant of the exterior medium. 
In the latter case, we find that for opposite charges, AF > always, implying a weakened attraction when compared 
to the baseline. For identical charges, however, there are cases for which the repulsion is actually enhanced compared 
to this baseline. Since it is possible to let e c — ► 1, this situation corresponds to an effective repulsion that is stronger 
than the vacuum case, a counter-intuitive result. We refer to this behavior as 'asymmetric screening'. 

When both repulsion and attraction are weakened compared to the e c baseline, which one is reduced more? This 
question is easily answered by considering 



5F = AF att - (— AF rcp ) = AF a < 



AF Ie 



When 5F > 0, there is a larger reduction of the attraction than of the repulsion, and vice versa. Using Eq. (|13p for 
AF rep and Eq. $T§§ for A.F att , we find 



SF = AF att + AF rcp = 2_ / fc e - 2fe (s+d) 



l~a 2 



1 — a 



L(i 



ae' 



-2kd\2 



? (1 + 



ae' 



-2kd\2 



dk > 0. 



For the case of the ej baseline, we see that AF is always negative for opposite charges. This indicates an enhanced 
attraction compared to the baseline (when both charges are in a uniform medium of dielectric constant £j). For 
identical charges we have AF > 0, implying that the repulsion is always enhanced when compared to this baseline. 
One can consider 

SF = AF att - (— AF rcp ) = AF att + AF rep . 



When SF > 0, the repulsion of identical charges is enhanced more then the attraction of opposite charges is. Using 
Eq. (HU) for AF rep and Eq. (21} for AF att , we find 



SF = AF att + AF rcp = 'L / ke- 2k ^ 
£i Jo 



(l + a) 2 



{l + af 



L(i 



ae 



-2kd\2 



' d ) 2 (1 + ae- 2kd ) 2 



dk > 0. 
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According to Fig. [5j asymmetric screening is quite pronounced at short ranges, and we expect the phenomenon 
to play an important role in biomolecular recognition and in the adoption of the native conformation of proteins. 
Particularly pronounced is the enhanced repulsion between charges of the same sign. This behavior should exert 
a rather strong veto on poor matching of charges as one part of a molecule interacts with another part or as two 
molecules interact with each other. Therefore, accurate calculation of electrostatic interaction is essential when 
considering biomolecular systems. 
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APPENDIX: SURFACE CHARGE METHOD 



The surface charge method[3]-[§| provides a relatively easy path to the induced surface charge. In the case of 
two identical charges, symmetry implies that the induced surface charge densities on the two surfaces are identical 
functions in the plane. Therefore we may write 

e \r- (d + s)z\ e c \r + (d + s)z\ J x , =+d \r-f>\ J z >^ d \r-r'\ 

The induced surface charge density a(p) is unknown, but can be expanded in a complete set of functions. Because of 
the cylindrical symmetry, Bessel functions are the obvious choice in this case. Any reasonably well-behaved function 
f(p) gives rise to the pair of transforms Q 



/>oo 

f(p) = / a{p)J v {fip)&0 
Jo 

poo 

a(P) = P / f{p)M0P)pdp, 



allowing us to write the surface charge as 

POO 

a{p)= / S(p)J v (pp)d0. 
Jo 

Furthermore, the denominator of the integrals in Eq. (|A.1[) can also be expanded in Bessel functions [lj|: 

V / dke im ^-^J m (kp)J m (kp')e-^ z >- z <\ 



m— — oo 



where z> = max{z, z'} and z > = min{z, z'}. 

In the vicinity of the surfaces, the potentials of the point charges are 

— p — JJ-, rrr = — / dkJ (kp)e~ 
e e |r - [a + s)z\ e c Jq 

and 

_—_——— = —/ dkJ Q (kp)e-' 
e c \r + (d + s)z\ e c J Q 
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The potential near the boundary at z = d due to the induced surface charge at z = d is 



z'=+d 



°° p CO 

V / dke m ^'^J m (kp)J m (k P ')e' k ^>- z ^ 
/ d/3S((3) / dfce- fc ^>- 2 <) V J m (fcp) / 
dp'p'J u (/3p')J m (kp') 



1(4,-4,') 



2ir / dpS(P) / dfce - *(*>-*<'j (fcp) 
Jo Jo 



dp'p'J v {f3p')Jv{kp') 



Letting v = turns the // integral into a standard one, [If 



and therefore 



Mf3p)MP'p)pdp = S MP) ( v > _i/2), 



-^^-dS' = 27r / dke- k(z >~ z <*>Jo(kp)S(k)/k. 



z'=+d \l " 

So for z' = +d and z > d (just above the top interface) 

poo 

2tt / dfc e - k{z - d) J {kp)S(k) Ik. 
Jo 

For z' = +d and z < d (just below the top interface) 

/>oo 

2tt / dke- k( - d - z) J {kp)S(k)/k. 



For z' = — d and z near 0? one finds a similar formula that is valid either above or below interface: 



2tt / dke- k(z+d) J (kp)S(k)/k. 
Jo 



The boundary condition at z = d is 



<9$ 



JS<G 



= e R 



z—d 



dz 



1 9z 

for every value of p, which leads to an equation easily solved for S(k): 

qkae- ks (e- 2kd - 1) 



S(k) 



2-Ke e (l - ae- 2kd ) 



Therefore 



o-(p) = / Jo{kp)S(k)dk 
Jo 

f °° /, s A;e- fcs (e- 2W -l) 
J (fcp)- 



27T£ C Jo 

qa 



(1 - ae- 2fed ) 



d/c 



2tt£ c J q 

qa 
2ire c 

qa 



, , (a-l)e~ 2kds 
J (kp)ke- ks 1 + ¥ ] dfc 



2ire r 



(1 - ae~ 2M ) 

/> DO 00 

J (M^e" fcs d/c+ / J (fcp)fce" fc(s+2<i) («- 1)^ a 11 e - 2knd dk 

r CO °° /*CO 

/ Jo(M^ fcs d/c + ^V l (a- 1) / J Q {kp)ke~ k(s+2{n+1)d) dk 
Jo n=0 Jo 
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Since all variables are real, we make use of the following integral [l6| 

/>oo 

/ e~ ax J v {j3x)x u+1 dx 
Jo 



„ +1 Aat _ (2a)(2/3rr(^ + (3/2)) 

V^(a 2 +/J2)^+(3/2) 



for v > -1 and a > 0. Recall that T(n + (1/2)) = ^(2n - l)!!2 _7l 3 so that T(3/2) = V^/ 2 - Therefore, the surface 
charge density is 



qa 
2ne c 



( s 2 + p 2)3/2 



X»-i) 



s + 2(n + l)d 



n=0 



[{s + 2(n+l)d) 2 +p'< 



i3/2 



from which it is easy to verify that J a(p)2npdp = 0. 

This charge density can be used to recover same energy and force as before. To compute the energy (and then the 
force), $(p = 0, z = s + d) must be computed from a(p). For z' = d, p = 0, and 2 = s + d, \f — r* 
Therefore 



*|2 _ s 2 +p ,2 



dS" = 2tt 



/" ^) dp- 



The p 1 integral is found in tables [17] to be 



xJ (xy) e ay 
■ da; = . 



(o 2 + a a ) 1 /a 



and so 



-e Jo 



-2fcs 



,-2fc<Z 



1 — ae 



-2A-r/ 



dfc. 



For z' = — d, p = 0, and z = s + d, \r — r*| 2 = (s + 2d) 2 + p' 2 . The contribution to the potential from the induced 
surface charge at z' — — d is 



^ldS' = ^ 



-e JQ 



-2k(s+d) _ 



-2kd 



1 — ae 



-2kd 



dk. 



The potential at p = and z = s + d is 



$(p = 0, 2 = s + d) 



-fc(2s+2<£) 



-2ks^ e -2kd _ ]\ a( ,-k(2s+2d) ^ £ -2kd _ j\ 



1 — ae 



-2kd 



1 — ae 



-2kd 



dk 



-o Jo 



2k d 



V_ / -2k(s+d) 1 ae 



1 — ae 



-2kd 



dk, 



and therefore 



-c Jo 



-2k(s+d) 



1 — ae 



2kd 



1 — ae 



-2kd 



dk, 



in agreement with Section ITTT1 Because U agrees, everything that follows from U must also agree. 
For opposite charges 



4> 



d \r-r> 



e c \r — (d + s)z\ e c \f + (d + s)z\ 
However, by symmetry er + = — er_ = a. The boundary condition yields 

qkae- ks (e- 2kd + 1) 



dS' 



^4d5' 



S(k) 



27re e (l + ae- 2kd ) 
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The surface charge density becomes 



s + 2(n + i)d 



<?{P) 



qa s 



27T£ C ( S 2 + p 2 ) 3 / 2 



[(s + 2(n+l)d) 2 +p' 



,2l3/2 



Again, it is easy to verify that J a(p)2npdp = and that the energy U reproduces the result in Section [TV] 
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